Diffusion, peer pressure and tailed distributions 
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We present a general, physically motivated non-linear and non-local advection equation in which 
the diffusion of interacting random walkers competes with a local drift arising from a kind of peer 
pressure. We show, using a mapping to an integrable dynamical system, that on varying a parameter, 
the steady state behaviour undergoes a transition from the standard diffusive behavior to a localized 
stationary state characterized by a tailed distribution. Finally, we show that recent empirical laws 
on economic growth can be explained as a collective phenomenon due to peer pressure interaction. 
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Fluctuations - measured by deviations from the mean 
value of an observable - of large systems often satisfy a 
Gaussian distribution. A classic example is a linear dif- 
fusion process Jl|-Q] which has numerous applications in 
many branches of science ||. However, there are many 
situations in Nature in which the probability of occurence 
of a fluctuation of size | A| is proportional to exp(— |A| P ), 
with p taking on a value greater than or equal to 1 (ex- 
ponential tails), but smaller than 2 Examples in- 
clude the temperature distribution of a Rayleigh-Benard 
system disordered systems such as foams || and 

glasses or granular materials [OJ, with even fatter tails 
in financial data 10,01. 



We present here a physically motivated advection 
equation and its exact steady state solution. The equa- 
tion has a drift term, originating from a kind of a peer 
pressure, of the same nature as that due to mechanisms 
of chemotactic signalling by microorganisms ||l3,[l4| or 
by the onset of cooperation in social groups [|15]l(||, or 
by competition between economic units |r7fl . We show 
that, on varying a parameter, there is a transition from 
diffusive behavior to a localized stationary state charac- 
terized by an exponential distribution. 

Consider a diffusional process in which random walk- 
ers move either to the right or to the left randomly and 
with no bias. In the long time limit, the distribution 
of walkers becomes flat and infinitely spread out. The 
mechanism for collective self-organization in our model 
is a kind of peer pressure. The non-linearity arises from 
an interaction (of spatial range £) between the walkers, 
which leads to a drift or a bias term which opposes the 
diffusional spreading and promotes aggregation. The ba- 
sic idea is one in which a walker perceives the populations 
of other walkers over a range £ both right and left of her 
own location and has a drift in the more crowded direc- 
tion. In the limit of small £, one obtains regular diffusive 
behavior, whereas a new class of steady state behavior 
characterized by non-Gaussian distributions is obtained 
when £ is sufficiently large. 

This idea may be encapsulated in a nonlinear equation 
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where P(x, t) is the distribution function of the locations 
of the random walkers at time t with a drift velocity 
v(x,t). The nonlinearity arises because v is itself a non- 
local function of P and is given by 



where 



v(x, t) = 
*+(x,t) = 



2 $+(x,t) +$-(x,t) 
dye {y - x) ^P{y,t) , 
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and A sets the scale for the drift velocity. Physically, 
$+(x,f) and <$>—(x,t) are measures of the population, 
within a range £, to the right and left of location x and 
the drift velocity is then a normalized imbalance between 
these two populations. 

In colonics of social individuals, each member can re- 
lease chemical substances in the environment (for exam- 
ple, pheromones) revealing the presence of food sources 
in a given area. The individuals in the colony are able to 
detect higher concentrations of the chemical signal and 
are attracted toward the food and move accordingly. The 
plausible assumption that the concentration of the chem- 
icals in a given region is proportional to the local den- 
sity of individuals translates into an effective interaction 
between members of the colony, and leads to the non- 
local drift term in Eq. (|l|). The interaction causes the 
net migration of an individual in the direction of the 
higher local density. The length scale, £, is a measure 
of the range of the biological sensory system of individ- 
uals. Thus Eq. ([l]) describes Brownian motion with a 
bias that mimicks attractive chemotactic signalling |l3| ] 
and promotes aggregation leading to a drift controlled by 
the population difference in localized regions on the right 
and left of x. Physically, the standard Gaussian solution 
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is obtained when there is no drift (A = 0) or when the 
interaction length scale goes to zero (£ — > 0). 

The existence of a transition in the steady state be- 
havior is revealed by a linear stability analysis of the 
homogeneous state. One can look for solutions of Eq. (|l|) 
in the form P ~ exp[ikx — j(k)t] representing spatially 
periodic perturbations. A direct substitution into Eq. (Fy) 
provides the exponential growth rate 



where 



7(fc) = D(l - r)k 2 
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Hence the homogeneous state is stable ("f(k) > 0) when 
r < 1. For r > 1 on the contrary, small perturbations 
grow with time and the homogeneous state becomes un- 
stable. Thus r plays the role of a control parameter of a 
phase transition. 

A mapping to an Hamiltonian dynamical systems al- 
lows us to derive the stationary soltuion of Eq. ([l]). In- 
deed the steady state distribution (P s = 0) satisfies the 
ordinary differential equation 
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where the prime denotes the derivative d/dx. The steady 
state values of $ + and in turn, satisfy 



= - p s 
= -<$>-/i + p s 
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obtained by taking derivatives of Eqs. (g) and (Q) re- 
spectively. It is convenient to introduce the new variables 
Q = + andn =$+-$_, so that Eqs. (§) and (§) 
may be written as 



n' = Q/^-2P s . 

Then, from Eqs. (g) we have 

Ek = AH 

P, 2DQ 
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which, with the aid of Eq. (|l(j) can be integrated with 
respect to x, yielding 



In P s — r In Q — const. 



(13) 



Thus, the steady distribution is P s (x) = CQ(x) r , where 
C is a normalization factor. 

The system of ordinary differential equations (Eqns. 
|i"o| and [ll]) can be regarded as an integrable dynamical 
system (whose energy is conserved) 



Q' = n/C 

IT = Q/e - 2CQ r 



(14) 
(15) 



describing the motion of a particle in one-dimensional 
potential 



Q 2 , 2C ^ r+ i 
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The spatial coordinate in the original system, x, plays 
the role of time in the dynamical system. The equation 
of motion can be derived from the Hamiltonian, 



H(U,Q) = ^ + V(Q) 



(17) 



which is a constant of motion, H(H, Q) = E in the dy- 
namical system. Its constant value E is fixed by the 
boundary conditions on P s (x). For instance, the require- 
ment that P s be normalized for infinite systems, implies 
that Q(x = ±oo ) = and H(x = ±oo) = 0, so that E is 
automatically set to zero. 
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FIG. 1. Plot of the potential V(Q) for A = 10, D = 1/4 
and £ = 1/31 (r = 20/31 - dashed line) and A = 30 D=l/4 
and £ = 1/20 (r = 3 - full line). 

This mechanical analogy allows one to obtain a com- 
plete understanding of the model with the key factor be- 
ing the shape of the potential energy (Eqn. |l^). A 
Qualitative change in the nature of V(Q) occurs (see Fig. 
[j]) when r crosses the value 1, i.e. for A£ = 2D, signalling 
a transition in the system's behaviour. When r < 1 
(dashed curve of Fig. |l|), only the solution Q(x) = 
is physically acceptable and thence P s (x) = 0. This cor- 
responds to a probability distribution which spreads out 
and vanishes as in standard diffusion. When r > 1 a 
non-trivial trajectory of the Hamiltonian system is pos- 
sible (full line of Fig. |l|) with Q vanishing when the fic- 
titious time x — > ±oo. This trajectory corresponds to a 
non-trivial asymptotic distribution (see Fig. g). The so- 
lution for r > 1 is obtained by integrating the differential 
equation 



Q' = i/CVQ 2 -gQ r+ \ 



(18) 
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which is obtained by solving for II as a function of Q (for 
E = 0) and substituting the result into Eq. ([Til), with 
g = 4C£/(1 + r). This leads to 



P s (x) = A{r) cosh 



L 2£ 



-(a; - xq) 



-2r/(r-l) 



(19) 



where A(r) = (r-l)T[(3r-l)/(2r-2)}/{2^T[r/(2r- 
2)]} enforces normalization and Xq is a constant which 
arises from the integration of Eq. ( |l8| ) . 

The presence of a single control parameter in the model 
can be deduced by a dimensional argument. Indeed, 
through the rescaling: X = x/£ and r = Dt/£ 2 , the 
coefficient A/2 in Eq. (||) becomes r and the dependence 
on £ disappears. Indeed when expressed in terms of X, 
the stationary state distribution depends only on r. 

The change of behavior ar r = 1 is a proper phase tran- 
sition. Indeed it is accompanied by spontaneous breaking 
of the translation symmetry x — > x + a of Eq. (|l|) for 
r > 1. The location, x$, of the distribution's center is 
dynamically selected and it depends, in a complex man- 
ner, on the initial conditions. Furthermore, one may de- 
fine a localization length, £, which is the spread of the 
distribution P s (x), and which diverges as I = |1 — r\~ v 
with a critical exponent, v = 1, signalling the transition 
to the dclocalized, translationally invariant steady state. 
The dynamics of relaxation to the steady state is char- 
acterized by an exponential decay with a characteristic 
time diverging as one approaches the transition. 

Another interesting feature of the system is the pos- 
sibility of spatially periodic steady-state solutions of 
Eq. (Q) for suitable boundary conditions, corresponding 
to periodic orbits of the associated Hamiltonian systems 
with negative energy. However, unlike Eq. (19), such 



solutions are expected to be unstable to perturbations. 

Many examples in which aggregation processes com- 
pete with diffusional tendencies can also be found in 
economics. Competition between economic units may 
lead to effective "peer" pressure captured by the drift in 
our model, whereas the diffusion term accounts for in- 
dividual idiosyncratic behavior. Such a situation could 
arise in charitable giving, when it is not anonymous. A 
given individual would, in an attempt to keep up with 
the Joneses, adjust his contribution in the direction of 
higher or lower giving depending on what his peers are 
doing. Of course, in this instance, the range of interac- 
tions, £, is limited because it is not feasible for an ordi- 
nary individual to mimic the behavior of a Bill Gates. 
In other contexts, tent-shaped distributions, which devi- 
ate sharply from a Gaussian, have been observed for the 
growth rates of firms, nations' gross domestic product 
(GDP) and complex organizations ( jlTj and references 
therein) . We have analyzed the time series of the growth 
rates of commodity prices across different Italian cities 
(see Fig. ||) and we find similar exponential tails. 
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FIG. 2. The scaled asymptotic distribution P s for A = 1 
and several values of r. Solid lines refer to the exact solution 
Eq. (|l!]). The open circles are obtained on numerically in- 
tegrating the differential equation (jjj) in the large t limit for 
r = 3. 
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FIG. 3. Histograms of the monthly growth rate R of 
two consumer price indices collected among 78 Italian 
Provinces within the period January 1995 - June 2001. 
The data are from "Italy's National Statistical Insti- 
tute" (ISTAT). The filled circles refer to expenditures for 
food and non-alcholic beverages, whereas the open circles 
show the data for miscellaneous goods and services (see 



tittp: //www. istat . it/Anotizie/Acom/precon/ indiceistat/paniere . ht 



for the entries). Both sets of data show the typical 
"tent-shaped" distribution in sharp contrast with the pre- 
diction of a Gaussian shape of Gibrat theory. The full and 
dot-dashed lines indicate a fit with Equation (|li^). 



The growth rate of an economic indicator b - such as 
the turnover of an organization, a nation's GDP, or the 
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price of a commodity in a city - is denned as 



R(t) = In 



6(* + 1) 
b(t) 



where t and t + 1 represent two consecutive time units 
(days, months, years etc.). The classical theory, due to 
Gibrat predicts that R(t) has a Gaussian distribu- 
tion. Empirical data, however, shows a tent-shaped dis- 
tribution on a logarithmic scale. In the case of firms, 
Stanley et al. jT^,[l9| have proposed an explanation in 
terms of the internal hierarchical organization of non- 
interacting firms. 

Our model provides a different explaination in which 
it is the interaction among firms which is responsible for 
the deviation from a Gaussian distribution. In a way, 
the peer pressure mechanism can be regarded as a re- 
gression towards the mean effect pp[ |. This term usually 
refers to an explicit "attraction" of R(t) towards a slowly 
varying moving average. Such an effect, however, is un- 
able to reproduce tent shaped distributions without ad 
hoc assumptions on the precise nature of the attraction 
potential. 

Here, instead, we refer to a collective phenomenon of 
"regression towards the population mean" : Each unit is 
attracted towards the instantaneous population average. 
It is the effective interaction among agents that generates 
a sort of "feedback" that induces the regression towards 
the mean mechanism. From a microeconomic point of 
view, it is reasonable to assume that the behavior of an 
economic agent is influenced by other agents. For ex- 
ample, a firm will strive to increase its growth rate if 
other peer firms grow at a higher rate. This interaction 
introduces correlations in R(t) that destroy the diffusive 
features and produces the characteristic exponential tails. 

In summary, we have introduced and studied, both 
analytically and numerically, a one-dimensional diffu- 
sion equation with nonlinear and nonlocal features. The 
asymptotic evolution of the equation leads to the walk- 
ers being attracted toward a state which exhibits a well 
defined non-Gaussian distribution characterized by ex- 
ponential tails independent of the nature of the initial 
condition. This mechanism of peer pressure may provide 
the basis for the development of more realistic models 
of self-aggregation and self-organization in cooperative 
states of populations of interacting individuals [Ell . 
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